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Localized modes of large amplitudes in nonlinear lattices are considered. 
The applied method allows the reduction of nonlinear problem to a linear 
inverse problem of phonons scattering on a local potential. The method is 
efficient in the case of strongly localized modes. Analytical description of 
such modes in monatomic chain is given. Results of numerical calculations 
of anharmonic local vibrations of light ions in pure and impure alkali halide 
crystal are presented. It is found that, in the case of amplitudes >0.5A, the 
vibrations depend very strongly on the crystallographic directions. 
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I. INTRODUCTION 

Recent studies of nonlinear vibrational dynamics lead to observation of long-living local- 
ized vibrations in perfect anharmonic latticeJ^ ( see also revie and references therein). 
The methods, which have been used so far for studying of these novel excitations, are based 
on direct numerical integration of classical equations of motions, being here the nonlinear 
differential equations. Due to the fast growth of the number of numerical operations and 
computational time with increase of the number degrees of freedom, the research is mainly 
focused on simple models, such as one-dimensional one- and two-atom lattices; only little 
have been done for real 3D crystals!. On the contrary, in harmonic approximation calcu- 



lations can be performed also for macroscopically large pure and locally distorted crystals. 
Here we present a method which makes use of harmonic approximation results for calcula- 
tions of the anharmonic local modes (ALM). The developed techniques allow one to reduce 
the nonlinear problem of ALMs to a properly formulated problem of linear local dynamics. 
It is applicable for lattices of arbitrary dimension, being especially efficient in the case of 
strongly localized modes, when it allows to obtain analytical solutions. The consideration 
is based on classic mechanics; account of quantum effects is given inS'&^l. 

The starting points for the method are the stability conditions of an ALM with respect 
to small fluctuations: 

• all frequencies of the spectrum of these fluctuations should be real and positive, 

• spectrum of fluctuations should contain also the basic frequency uJi of the ALM, i.e. 
there should be harmonic local mode (HLM) with the same frequency uji\ the same 
holds for the higher harmonics of the mode, 

• the ratios of the amplitudes of atoms of the HLM should coincide with that of corre- 
sponding ALM. 

The physical meaning of the condition a) is obvious: the conflguration of the crystal 
lattice with the ALM must be stable (not metastable). The conditions b) and c) ensure 
that a solution, which describes an ALM are stable with respect to inflnitesimal changes of 
the ALM phase. The last two conditions allow the self-consistent calculation of the ALMs. 
Indeed local mode of flnite amplitude in anharmonic lattice causes local perturbation of the 
dynamical matrix of small vibrations. Therefore, knowing (or choosing as a trial parameters) 
the amplitudes of the ALM on different atoms, we can determine the changes of the elastic 
springs, caused by the ALM. Then, by applying the Lifschitz formula of local dynamics in 
harmonic approximation, the expressions for the HLM (and ALM) may be obtained, both for 
its frequency and amplitudes of contributing atoms. These expressions allow the obtaining 
of the self-consistency conditions, being the set of of equations for the parameters of the 
ALM. Thus, the method reduces the nonlinear problem of calculation of an ALM to a kind 
of a linear inverse problem of phonons scattering on a local scattering potential for small 



vibrations (phonons). If the ALM is strongly localized, then the procedure is simple while 
one needs to calculate only few amplitudes. This can be done analytically. The method is 
applicable for lattices of different dimensions. Besides, it can be also used for calculations 
of anharmonic local modes of defects. 

To illustrate, how the method works, we will apply it first to the ALMs in monatomic 
chain with hard quartic anharmonicity. We will show that one gets analytical results which, 
in the limiting case of high frequency, coincide with known solutionJl. Then we will calculate 
odd ALMs in pure and impure alkali halide crystals with light and heavy ions and with 
account of quartic anharmonicity (the ALMs under consideration are almost fully localized 
on the light atoms). Frequency spectra in these crystals are calculated in the shell mo 
Anharmonic constants are determined from the Born-Coulomb-Mayer potentials. 



II. GENERAL 

Let us consider a vibrational system with the potential energy 

^ nin2 nin2n.3 

where f/„ are the Cartesian displacements of atoms situated at the site n, V^'^^ , V^^^ , . . . are 
harmonic and anharmonic springs. The displacements [/„ satisfy the following equations of 
motion 

- M„f/„ = y: v;Sf/„, + E ySUUn,u^2 + • • • , (2) 

ni n\n2 

where M„ are the masses of atoms. We suppose that a ALM is excited at the site n = 
and its nearest neighbors. This excitation is the solution of the equations of motion with 
exponentially decreasing \Un\ with \n\ ^ og. As it is knowniSB a solution, which corresponds 
to the ALM with the frequency ui has the form Unit) = An cos{uJit)+^n + 0{uJi), where 0{u!i) 
is the sum of small terms of the frequencies 3uJi, 5uji, . . . (i.e. the higher order harmonics; 
these harmonics are neglected below). The amplitudes An and the shifts ^„ can be found 
by equating the coefficients before the terms with the same time dependence&SQ. The 
corresponding equations are 



nl n\n2 

+ 3 E Kinln2n3 (7^ni^n2 + ^nxin^^Un-j, + • • • , (3) 
nnin2ri3 

E KiSe-i + E K!nln2 (|f^nif^n2 + UU) + • • • = 0. 

ni ni7i2 

To take into account small fluctuations we present the displacements in the form Un = 
Un + (Inl \/^ni where g„ stands for the reduced small displacement. The equations for g„ 
are obtained from and read: 

'in = E (^2nn' + W2nn')qn', (4) 
n' 

where V2nn' is the dynamical matrix in harmonic approximation, 



+ ^ E V^X^.iUn^Un, + 2UU) + •••)' (5) 

n\n2 

describes anharmonic renormalization of the dynamic matrix due to the ALM (time- 
dependent terms ~ cos kuit, k = 0, 1,2, . . . are neglected here; these terms are important 
in quantum mechanical description of fluctuations: they lead to multi-phonon decay of the 
ALMii). In the case of strongly localized modes one can include into consideration only 
small number of amplitudes Then {VF2rm'} is the A^xA^-matrix of the local springs of 
the lattice. 

Hamiltonian, which corresponds to the equations of motions (^, equals Hp^ = -ffo,p/i + ^2, 
where 

Ho,ph = ^ E + ^ E ^2nn'gngn' = ^ E(^? + i^) 
nn' i 

is the phonon Hamiltonian of the harmonic lattice, Xj = J2n ^inQn are the normal coordinates 
of the lattice (in a perfect lattice ~ e**^"), 

^2 = ^ E ^2nn'qnqn' 
nn' 

describes the effect of the ALM on the phonons, W2nn' are corresponding local distortions 
of the elastic springs, given by (Bl). 



The HLMs, induced by the interaction V2, are determined by the poles of the spectral 
Greens functions 

POO 

Jo 

on the real axis of iv; Gnn'{t) = Si (^in(iin' sin(a;jt) is the retarded matrix-Green's function 
of phonons, ei„ = Cin/cOi . This matrix- function can be found by applying Lifshitz formula 

Gnn' = G^nl' + X] G^nli^2nin2G n2n' , (7) 
nin2 

where G^l,{uj) are the Green's functions of the harmonic lattice. These functions can be 
calculated by standard methods of lattice dynamics. To find the frequency and shape of the 
HLM (and ALM) let us introduce the configurational coordinates 

n V 

The first transformation S is chosen to diagonalize the perturbational quadratic form V^'- 

V 

The second transformation s{uj) diagonalizes the dimensionless Green's function matrix 
Gyy' = y/r]yr]i,'G'^!!li'- 

G\iIi = ^mm' y^.S^vS^i,iGvi,i = 6fj,n' Rp,nRp.n'Gll,. 
vv' nn' 

Here G^^li{uj) are the lattice Greens functions in Qi^-space, i?^„(ci;) = Y,v Snv{^)y/>li'Sun- In 
the (5//- representation the Lifshitz formula takes the form G^^{uj) — Gf^^{ijj){l — G ^^{u;)) . 

Prom symmetry considerations one can chose the configurational coordinate Q/i which 
corresponds to the ALM under investigation. The frequency uji of the local mode is given 
by the position of the pole of the Green's function G^^(a;) on the real axis cu, i.e. satisfies 
the relation 

= 1; (8) 

the relative amplitude of the n' atom in the mode is given by the formula: a^n = 

-R/in('-^0'\/^/I^A'M where G^^^' (cu) = dG^^^iu) / duo . Taking into account that the rel- 

ative amplitudes satisfy the normalization condition J2n(^l.n 1, we obtain: 



The relations (y) give the set of equations for a^„. In combination with (8) they allow one 
to find the amplitudes An of the ALM of given frequency ui. The method works better for 
strongly localized ALMs with high frequency, when the number of contributing amplitudes 
An (and equations to be solved) is small. 

In the case of high frequency of the ALM (tu/ ^ ujm, ujm is maximal phonon frequency) 
the non-diagonal elements of the Gnn'i^i) are much smaller than the diagonal ones. Then 
G^^^,{uJi) ~ 6i,,yiGQQ (uji), i.e. the Greens-function matrix diagonalizes simultaneously with 
V2. In this case the normalized amplitudes of the ALM can be determined by simple self- 
consistency condition = Syn', absolute amplitudes one can find from the relation Guu{uji) = 

If the limited number of atomic displacements is taken into consideration then the 
quadratic form V2 does not depend on the totally-symmetric coordinate, being the com- 
mon displacement (the sum of equal displacements of all accounted atoms). This means 
that all modes, which contribute to the diagonalized quadratic form V2, are orthogonal 
to this displacement. Evidently, the same holds for the liner combination of shifts, which 
give contributions to the ALM. Otherwise, the sum of displacements of all atoms, which 
contribute to the ALM, should be equal zero: 

E^n = 0. (10) 

n 

The same holds for the sum of momenta of atoms at any time- moment. This property of 
the ALMs may be called as immobility condition. 

III. EVEN AND ODD ALMS IN MONATOMIC CHAIN 

We apply the method to description of ALMs in monatomic chain. The potential energy 
of the chain with account of the nearest-neighbor interactions has the form 

!/('■) 



In harmonic approximation (K^^^ > 0, V^"^^ = 0, r > 3) for uj/um > 1 the Green's 
functions of the chain equafl 



^(0) ( ^ _ (-P)'" "''^Af/W 



2 
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ojm = 2-^/^(2) /Mo is the top phonon frequency, p = {uj/um — \J^'^ /^"m — 1)^ < 1- In the 
limit of strong amplitudes and high frequency one gets two strongly localized ALMsi'i: 

1. the odd mode with An = A_n, Ai ^ —Aq/2, \An\ <^ \Ai\, n> 1 and 

2. the even mode with An+i = —A_n, |^n+i| |Ao|, n > 1. 

For 72 ^ 1 we have A|n|+i — —pA\n\- For smaller frequency the wings (tails) of the mode are 
stronger. Below we give analytical description of strongly localized ALMs, having remarkable 
uji/ujm — 1, in the chain with quartic anharmonicity (r = 4). Note that in the opposite limit 
of very small uji/ujm — 1 there exists an analytical solution of the problem, which describes 
the ALM of large size!. 



1. Even mode 

Let us consider first the even mode. With account of four central atoms displacements 
the contribution of the ALM to the potential energy of phonons has the form 

V, = K,Al[b{{q, - q,Y + (go - q-i?) + (gi - go)'], 

where K/^ = QV^'^yMQ, 6 = (1 + /5)^/4, /3 = A2/AQ. Two following even modes contribute to 
the V2 and ALM: 

yi = qo-qiV^, y2 = q2-q-iV2. 

In the strong localization limit the main contribution is given by yi. Contribution of y2 
depends on the ui being larger for smaller ui. Our task is to describe the last dependence. 
In yi, y2— subspace V2 = K4{{2 + b)yf — 2byiy2 + by"^)- This quadratic form is diagonalized 
in the rotated basis Qi = yi cos (p + y2 sin 0, Q2 = —yi sin 4> + y2 cos with tan (20) = b. 
In the limit of high frequency uji the self-consistency condition reads 

sm = /? ^ - = ^ —. 

^ ^ 2 8 



One gets /5 1/6 in agreement with the corresponding calculation of (3 on the basis of the 
nonlinear equations of motionll. To describe cu;— dependence of (3 we should find coordinates 
(5i,2 which diagonalize the Green's function matrix G^^^, = Gqq {y^rj^Aj^guiy'), = 1,2, 

gu = Pu cos^ (f) + p22 sin^ (j) + pi2 sin (20), 
922 = Pu sin^ + P22 cos^ - pi2 sin (20), 
5^12 = -pi2 COS (20) + ^(pii - P22) sin (20), 

where pn = 1 + p, P22 = 1 + P^, P12 = p(l + p)- The diagonalization is achieved by rotating 
of the coordinates (5i,2 on the angle a with 

tan (2a) = 7 — -. (11) 

{Vi9u - V2922) 

The self-consistency condition now is 

/3 = sin(a + 0). (12) 

This is equation for l3{uJi), which can be easily solved numerically or (approximately) ana- 
lytically. E.g. in the strong localization limit, when both, a and are small, one gets 

P^- + ^^ . (13) 

6 12v^(l + p) ^ ^ 

The dependence of the dimensionless frequency of the mode Qi = uji/um on the amplitude 

Ai is given by the equation (8) with u = 1: 



gS?(^0 = K.AlG^Si^i) (l + b + Vr^) \9ii cos^ a + 922- sin^ a + 912 J- sin {2a] 



= 1. 
(14) 



For uoi > 1.15c(jjv/ given formulas describe the ALM rather well: the contributing amplitudes 
of the next nearest atoms is less than 10~^, thereat their contribution to the energy of the 
ALM is less than 10"^. 



2. Odd mode 

Let us consider now the odd ALM with account of displacements of five atoms. The 
mode under consideration satisfies the parity and immobility conditions An = A^n, A2 + 



Ai + Aq + A_i + A_2 = 0, which give Ai = -Ao(l/2 + f3), (3 = A2/A0. The contribution of 
the ALM to the potential energy of phonons equals 



V2 = K^Al 



Kfe - qif + (g-2 - g-i)') + (gi - go)' + (g-i - go)' 



where = 9VW(1 + 2/?/3)V8, 6 = (1 + 4/?)V(3 + 2/?)^. Two odd modes, which give 
contributions to the ALM read: 

yi = :y^(2go - gi - g^i), 1/2 = y^(3g2 + 3g-2 - 2go - 2gi - 2g_i). 



In the ?/i^2-space V2 = K4[{3 + b/3)yf — {2^/5/3)byiy2 + (5/3)6y|]. This quadratic form is also 
diagonalized in the rotated basis with sin (0) ^ 6/2, where b = 2v^/(9 — 46). In the limit 
of large uji/ujm the main contribution to the ALM is given by the ?/i-mode; contribution 
of the 1/2-niode is given by the self-consistency condition (3 ~ 3 sin (0)/2-\/5 ~ 6/6. Taking 
6 1/9 one gets the value [5 = 1/54 in agreement withll. A more accurate approximation is 
p ^ 3/131. 

The w;— dependence of /3 and Aq- dependence of ui can be found in the same way as for 
the even mode. In this case 

r^i 2 = _^iL(3V5 + 156 + (3^5 + gb))(l + ±Vl^). 
2V5 + 46 

For small f3 one gets 

The factors g^,^r are also determined by the same expressions as in the even ALM case but 
with the following p^^r. 

. 4p p2 

^ 15 15 5 5 ' 

Pi2 = v^p(^ + | + ^). 

The self-consistency condition in this case is: 

(3 = 3sin(^ + ^^ ^^^^ 

2v 5 cos {(j) + a) — 2 sin (0 + a) 



where a is determined by (0); thereat expressions for 771^2 and for p^^i are given above, g„^i 
are also determined by the same formulae as in the even ALM case. In the case /3 -C 1 one 



gets sin (a) ~ P\j'r\2l^'r\\ where r/i ~ 3 — 6/3, r]2 ~ 56/3, p = p(5 — 3p)/(3 — p). This gives 

The dependence of the frequency of the mode is given by equation (6) If uji > 1.15ujm then 
the odd ALM is well localized [jS < 0.065), more than 99.98 percent of energy of the mode 
come from 4 central atoms. 

The calculated dependences of fli and /5 on amplitude for even and odd modes are plotted 
on Fie;.!. 



IV. ALMS IN ALKALI HALIDES 

We present here calculations of odd local modes associated with light impurity and host 
ions in alkali-halide crystals. In these cases the ALM is almost fully localized on the light 
ion and the problem reduces to calculation of the frequency of the mode in dependence on 
its amplitude. 

Within the approximation of the nearest neighbors interaction the potential operator has 
the form 

00 1 

a H m=l 

where a = ±x, ±y, ±z are the directions to the nearest neighbors, n = {nx,ny,nz) is the 
vector of the lattice sites, is the vector of the site nearest to n in a direction, V^"^^ = Vm 
is the m-th derivative of the pair potential between atoms (ions) n and at their mean 
distance i^ona, Rn^ = [{Ron^ + fanj'^ + " rlf^J^^"^ - Ron^ Is the operator of distance 
between the nearest neighbors in the a-direction, = q^fi — q^fi- , qp is the /5-component 
of the displacement vector q^ of the atom n, a, P = x, y, z, f|_ = f^^_ + f^^_ + r^^_ . By 
expanding V in the power series of displacement operators fc^-one gets 



where 



2 aria ^ 



+ ...], 



(18) 



= ViR^\ Vi = (1/2 - V^)Ro\ 



(19) 



V2, V'i and V4 make account of the central, while V2, V^, VI and V^ of the non-central 
forces. The potential considered does not take account of the covalent interaction which 
leads to the chemical bonding. This (covalent) interaction can, however, be easily included 
in calculations by introducing additional terms of the type V2, V3 and V4. 

As it is known, harmonic non-central springs are normally 5 to 10 times smaller than 
the central springs V2. Our calculations of V3 and V^ for alkali halides show that V^/V^ is 
even smaller than V2/V2 {V^ is twenty-thirty times smaller than ^3). The same holds also 
for quartic and higher order anharmonic terms: the higher order anharmonicity the smaller 
are corresponding non-central interactions as compared to central ones. Therefore, as a first 
step, only central forces may be accounted. One can show (see Appendix) that the main 
effect of the non-central anharmonic interactions to the static local dynamics consists in 
renormalization of central elastic constants: U'^V^ is replaced by t/^(V4 -|- 6I4). This allows 
one to improve the central-force approximation by replacing V^hy V4 — V4 + 6V^'. 

We calculate a strong local vibration of light impurity or host atom (ion) situated at 
the origin of our reference frame. In this case solutions of classical equations of motion, 
corresponding to the local mode, satisfy the conditions: 1^40 1 » \Afi\. This allows one to 
suppose that the mode is well localized on the atom at the site n = 0. Then in the approxi- 
mation of central forces 9 coordinates of atoms contribute to perturbation of the dynamical 
matrix: 3 coordinates of the central atom and 6 directed to this atom coordinates of the 
nearest neighbor atoms. We chose coordinates according to Fig. 4. In this representation 
the impurity induced change of the dynamical matrix V2 is 



V = 



^ ^ 



Wy 







/ . _ _ \ 

Pa Ta Ta 



-la 7« 



In harmonic approximation f3 and 7 do not depend on a: 7 = AV2 is the change of the 
magnitude of the central elastic constants due to defect atom, /3 = 27 + — M/Mq), lo 
is the frequency of the normal mode, M/Mq is the ratio of impurity and host atom masses. 
Quartic anharmonicity leads to amplitude dependent corrections of elastic constants and to 
their dependence of a: 



,/ M ^ 
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(20) 



Ac is the a's Cartesian component of the amplitude of the local mode. Perturbed Green's 
functions are 
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where 



Gn 



G^:i = E 



^nvk 



Up' — tU^-' 
vk 



([/-G(°)ug-^G(°))„„, 
E 



^(0) _ V- ei^^e2^fcCos(fc^c/) 



vk 



(21) 



vk vk vk 

are Green's functions of perfect lattice, d is the lattice constant, k is the wave vector of 
phonon, v is the phonon branch, cu^g is the frequency of phonon, e^^^^ is the projection of 
the polarization vector of phonon onto a;-component of the atom. We use for the Green's 
functions G^^-* their values calculated within shell modellll. 

The frequency of the ALM, which depends on the amplitude of the mode, is given by 
the position of the pole of the Greens function Gaool^) ^^^^ ^^^1 ^^i^ of ^- Energy of the 
local mode depends on the amphtude via dependence of the uJi and directly: 



\muj]A] + U^,{A\ + + A^) 



(22) 



{A} = Al^Al^Al\ 



In a perfect lattice the perturbation matrix V ~ V^A^. This perturbation can lead 
to appearance of the ALMs. In three dimensional crystals they appear only if the local 
perturbation is strong enough. This means that there is an minimal (critical) amplitude of 
the central atom Aq for appearance of an ALM''. 



We performed calculation of the ALMs and ILMs associated with light F~ and Na~ ions 
in different alkali halide crystals. Results of some of our calculations are presented on Figs. 
3 and 4. One sees that frequency of the ALM of light ion F~ in KCl rather strongly depends 
on amplitude of the mode and on the crystallografic direction of the vibration. Thereat 
the smaller is the distance to the nearest neighbor {nn) atom in direction of vibration the 
larger is effect of anharmonicity at the same amplitude. The reason for this correlation 
is rather obvious: central anharmonic forces, which dominate in the anharmonic repulsive 
interactions, are stronger for smaller nn distance. One also sees that the ILM on light 
host ion Na+ in Nal appear already at rather moderate amplitude ~ 0.4A. Thereat with 
increasing of amplitude the mode first appear in (100) direction, then in (110) direction 
and then in (111) direction. The reason for such directional dependence is given above: 
anharmonic interaction is stronger for smaller nn distance. Note that the decay time of 
the modes due to two-phonon emission, as it was shown inS also strongly depends on the 
crystallografic direction. 

In conclusion, the method of calculation of anharmonic local modes both, ILMs (in pure 
crystals) and ALMs (in doped crystals) is proposed. The method allows the reduction of the 
problem of nonlinear local dynamics to the linear inverse problem of phonons scattering on 
a local potential, caused by the local mode. New analytical description of the even and odd 
strongly localized ILMs in the monatomic chain is given. Results of numerical calculations of 
anharmonic local modes of light ions in pure and impure alkali halide crystals are presented. 
It is found that the ILM on a light host ion appears already at rather moderate amplitude 
(e.g ~ 0.4A in Nal). Both, ILMs and ALMs in alkali-halide crystals strongly depend on the 
crystallographic direction. 
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A. Appendix: Non-central forces 



Although non-central anharmonic interactions are remarkably weaker than central ones 
their number is larger. Therefore it is of interest to account them. These interactions 
switch-on into the perturbation V2 all 21 Cartesian coordinates of the central atom and 
its 6 nearest neighbors. It is convenient to chose these coordinates as follows (see Fig. 1): 

qi = xo, q2 = yo, qs = Zq, q^ = Xi^, = yi^, gg = ^1., qr = X-i^, qs = y-i„ qg = z_i^, 
qio = hy, qu = Xi^, qu = Zi^, q^ = y-i^, qu = x^^, ^15 = ^-i,, qw = ^1,, qn = Xi^, 
qi8 = yi^, qi9 = ^-i^, ?2o = x-i^, 521 = y-u', U = ^Juf+Uf+U^. in this representation 
the perturbation matrix of lattice dynamics equals 

/ _ _ _ \ 

i/Q -Ui -V2 -U3 
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-pj i>i 

-pj h' 



—vJO z/i" 



where — (31^ is 3 x 3 matrix, PjC") = (1, 1) x v^c") are 3x6 matrixes, i>i — I2 ^ 1^1 are 6 
X 6 matrixes, /„ is nxn-unit matrix, 7 and 7' are changes of central and non-central elastic 

springs due to the impure central ion, /3 = ci;^(l — M/Mq) + 27 + 47', 
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In harmonic approximation — vy — i/j//. Quartic anharmonicity causes amplitude- 
dependent renormalization of the elastic springs and leads to the following corrections of 
(additions to) given above matrixes u: 



1. (lOO)-direction: u'y = u[„ = 
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where 6i = C/^(Vi + 6^4)74, S2 = Lf^Vl/A, Ci = 25i, (2 = 4^2; energy correction equals 
E'i^U\V4 + 6V;)/6A; 



2. (llO)-direction: u[, — v[, 
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where(Ti = [/^(T^ + 7^40/8,(72 = V^V'J^.a = 4(72,^1 = 2(7i + 6(72,6 = 6(72; energy 
correction equals E[ = U^{Vi + 12VI)/12^- 



3. (lll)-direction: i/J/ = u[,. 





' <^ f 




( 

Ki 


«2 


\ 




/ 

«2 


«3 


"1 




/ 

«2 


«3 







f <!> f 




K2 







, ^2 = 


Ki 


K2 


K2 




K2 


«3 


^2 








\ '^^ 





y 









1^3 J 






H2 


^^2 



/ 

where ki = U^{Vi + Wl)/12, K2 = K3 = K2/2, (p = 2ki + Ak2, f = 8^3; energy 

correction equals E'l = U^iV^ + 18\/4')/192. 

Corrections ~ V^' are neglected; they are at least one order of magnitude smaller than small- 
est accounted correction ~ V^. As one sees the main effect of the non-central anharmonic 
interactions to the local dynamics consists in renormalization of central elastic constants: 
A^V^ is replaced by A\V4 + 6Vi). 
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Figure captions 



Figure 1. (5 (solid lines) and Ai (dot-dashed lines) vs. uji for even and odd modes. We used 
for calculations ujm — 1 and y^^^Mo = 1. 

Figure 2. Atomic displacements, which contribute to the perturbation of the dynamical 
matrix: a) account of central forces - 9 displacements; b) account of central and noncentral 
forces - 21 displacements. 

Figure 3. Frequency of the localized mode (wi/wm) vs. amplitude (in A). KCl : F, [100] 
- solid line, [110] - dashed line, [111] - dot-dashed line. 

Figure 4. Frequency of the localized mode (wi/wm) vs. amplitude (in A). Nal, [100] - 
solid line, [110] - dashed line, [111] - dot-dashed line. 
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